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We present a short review of how the effec- 
tive action formalism, well known in 
relativistic quantum field theory, can be 
used to discuss Bose-Einstein conden- 
sation of non -relativistic gases. This method 
lends itself very naturally to an inter- 
pretation of Bose-Einstein condensation in 
terms of symmetry breaking. It also 
allows for the definition of a very elegant 
regularization technique involving 
generalized ^-functions. We show how this 
method can be used to recover the well 
known results for the free boson gas, as 
well as the charged boson gas in a 
constant magnetic field. A general criterion 
for interpreting Bose-Einstein conden- 
sation in terms of a phase transition with 



symmetry breaking is given. Finally we 
present an analysis of Bose-Einstein 
condensation in a harmonic oscillator 
confining potential trap, and show how the 
results of this simple model are in excel- 
lent agreement with experiment. 



Key words: Bose-Einstein condensation; 
effective action; symmetry breaking; zeta 
functions. 



Accepted: May 15, 1996 



1. Introduction 

It is now well over seventy years since the phe- 
nomenon referred to as Bose-Einstein condensation 
(BEC) was first predicted [1,2]. For s system of non- 
relativistic spin-0 bosons in three spatial dimensions, a 
discussion of BEC is now part of any undergraduate 
course in statistical mechanics. Until recently the best 
experimental evidence that BEC could occur in a real 
physical system was liquid helium, as suggested origi- 
nally by London [3]. However although the behavior of 
liquid helium at low temperatures can be qualitatively 
described by the free boson gas model, the detailed 
behavior deviates substantially from this simple model. 
Physically this is of course because the effects of inter- 
actions which are neglected in the free boson gas model 
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are important in liquid helium. More recently it was 
suggested [4,5] that BEC could occur for excitons in 
certain types of non-metallic crystals (such as CuCl for 
example). There is now good evidence for this in a 
number of experiments [6] . 

The most exciting experimental evidence for BEC has 
come from the recent observations of very cold alkali 
gases. BEC has now been observed to occur in gases of 
rubidium [7], lithium [8], and sodium [9]. These sys- 
tems are very dilute and as a first approximation are well 
described by a boson gas model with no interactions 
among the atoms. The atoms are confined in a magnetic 
trap which can be modelled by a harmonic oscillator 
potential. We have recently discussed [10,11] how the 
occurrence of BEC is to be interpreted in such a system, 
and given the details of the harmonic oscillator potential 
trap, are able to calculate a characteristic temperature 
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which is in excellent agreement with the values found in 
the experiments. This work will be reviewed briefly in 
Sec. 6. 

The main purpose of our review is to show how the 
effective action formalism may be used in a straight- 
forward and natural way to discuss EEC. The general 
formalism is presented in Sec. 2. In Sec. 3 we introduce 
a particularly useful method for defining the effective 
action, and show how the usual thermodynamic poten- 
tial may be recovered. The interpretation of EEC as 
symmetry breaking is given in Sec. 4, and used to study 
two models: the free boson gas, and a charged boson 
gas in a constant, externally applied magnetic field. A 
general criterion for deciding if EEC will occur as a 
phase transition with symmetry breaking is presented in 
Sec. 5. In Sec. 6 EEC in a harmonic oscillator confining 
potential is discussed, and we describe how this model 
compares favourably with the experimental results on 
alkali gases. 

2. General Effective Action Formalism 

In this section we wish to discuss briefly the effective 
action approach to quantum field theory at finite tem- 
perature and density. We will then use this formalism to 
see how EEC may be understood in terms of symmetry 
breaking, since this interpretation arises in a very 
straightforward manner within the effective action 
framework. One advantage of adopting the effective 
action approach is that it may be applied to situations of 
great generality, such as curved configuration spaces of 
arbitrary dimension, spaces with boundaries or with 
complicated topologies, or situations where background 
gravitational or electromagnetic fields are present. 
Furthermore, as we will discuss in the next section, the 
effective action formalism allows for a very elegant 
mathematical regularization procedure to be used. 
Finally, the extension of the method from noninteracting 
to interacting gases may be performed in a systematic 
way. 

Our attention will be on a system described by 
a nonrelativistic Schrodinger field ^ with action func- 
tional 

5[^, ^^]= \dt j da-J^(^^^-^^^) 



-^\D^\'-u,(x)\^n. (1) 

Im J 

(We have adopted units for which ^= \.) Here S repre- 
sents the spatial configuration space. It can be any 
Riemannian manifold, with or without boundary; 
however, as our applications will be confined to flat 
Euclidean space ^^,S may be thought of as a finite box 
in ^^ with periodic boundary conditions imposed on 



the sides of the box. The infinite box limit will be 
understood, dcr^ represents the invariant volume element 
on S, which for ^ = ^^ is simply given by dcr^ = d^x. 
D is the dimension of the space which we keep arbitrary. 
Ui(x) represents any potential, which is assumed to be 
time independent. D = V-ieA is the gauge covariant 
derivative, with A the vector potential describing any 
background electromagnetic field which might be 
present. We adopt the gauge choice 



Ao = 0, V • A = . 



(2) 



It is possible to add a self-interaction term for the 
Schrodinger field to Eq. (1), or generalize in other ways 
by considering a number of different fields. The general 
formalism of this subsection does not depend in any 
significant way on the precise form of the action Eq. (1). 
In addition to the action describing the Schrodinger 
field, we must include the action functional for any 
background gravitational or electromagnetic fields 
which are present. In this review we will only consider 
the case of time independent background magnetic 
fields. The action for the magnetic field will be taken as 



H 



da^ 



- F F^J- T^ A- 



(3) 



Here A ' are the components of A , and Fij = V,A ~ V^A, 
is the field strength tensor, /ext represents the compo- 
nents of the current /ext responsible for setting up the 
background magnetic field. If we have D = 3, then the 
magnetic field vector B with components B' may be 
defined by Fij = SfjkB ^, where 6^^ is the antisymmetric 
Levi-Civita tensor. However if we keep the spatial di- 
mension D general, the magnetic field is not described 
by a vector and we must deal with the antisymmetric 
tensor F^ . 

The Schrodinger field action Eq. (1 ) is invariant under 
the local gauge transformation 



A(x)->A(x)-hVO(x), 



(4) 



(5) 



where ^ (x ) is an arbitrary function of the spatial coordi- 
nates x on ^ . Associated with any local gauge symme- 
try is a conserved current (defined via Noether's theo- 
rem). In the case of Eqs. (1), (4), (5), the conserved 
current is 



J=^ (V^^ ^-^'^V^)--A I ^P . (6) 

2m m ' ' 
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The conserved charge associated with this current is 



Q 



./d.. 



I "^1' 



(7) 



Apart from a factor of e, Eqs. (6) and (7) may be seen 
to be the probabihty current, and probabihty of wave 
mechanics. Provided that we restrict the background 
electromagnetic field to sufficiently small values to 
prohibit pair production, we will have a conserved parti- 
cle number 



N=Qle= \da,\^ 



!• 



(8) 



In the case of neutral systems, such as atoms confined 
by a magnetic trap, it is the number of particles which 
is conserved. We will therefore consider A^, rather than 
2, to be conserved in what follows. 

The thermodynamics is described by the grand parti- 
tion function ^. In order to incorporate the conserved 
particle number A^ it is customary to introduce a 
Lagrange multiplier ^ called the chemical potential. 
Then we may write 



tr exp (-I3(H -fxN)) 



define the effective action in a way we will describe 
later. These sources are also useful for obtaining the 
Green functions of the theory. The periodic boundary 
conditions imposed on the fields in the path integral 
ensure that the Green functions obey the usual boundary 
conditions for finite temperature field theory. (See Ref. 
[15] for example.) 

It is important to keep in mind what is being held 
fixed when the path integral Eq. (10) is performed. The 
path integral is computed with the temperature T, vol- 
ume and metric on X, chemical potential fju, background 
gauge field A, and Schwinger sources Js and /J all 
regarded as fixed. By performing Legendre transforma- 
tions on the appropriate variables, it is possible to obtain 
other functionals which hold different variables fixed. In 
place of the chemical potential fju, we wish to keep the 
particle number fixed. In place of the Schwinger 
sources Js and /J, it_ proves advantageous to introduce 
background fields ^ and ^^. Because we will only 
consider_ the transformation from (fju, Js Jl ) to 
(A^, ^, ^) we will not indicate the functional depen- 
dence on any other variables which the partition func- 
tion depends on. 

We will define 



(9) 



w[fjL,Jsjn=-in^[fjL,jsjn 



(12) 



Here H is the Hamiltonian operator for the theory Eq. 
(1), and A^is the number operator which is obtained from 
Eq. (1) by regarding ^ and ^^ as field operators. 

We will use the path integral method [12,13] to com- 
pute ^. To do this it is easiest if we adopt the imaginary 
time formalism, where the path integral extends over all 
fields which are periodic in imaginary time with period 
13= l/(kT). The grand partition function may be ex- 
pressed as 



-I 



^= [d^'^d^]c- 



(10) 



where 



S = 5ern+ | dt Ida, U (^"^4^- 4^^^) 



2m' I '^ I I II 



+ /Jl/^+ ^"^Js 



(11) 



is obtained from Eq. (1) by performing the Wick 
rotation t -^-it to imaginary time, and including the 
conserved particle number Eq. (8) with the Lagrange 
multiplier ^t . A Schwinger [ 1 4] source Js and its com- 
plex conjugate Jl have been introduced in order to 



In interacting quantum field theory W is the generator 
of connected Green functions. We will eliminate the 
dependence on the Schwinger sources Js and j] by 
defining 



'^— 



8W 
8/1 



^Js 



(13) 



8/. 



Ai,/J 



(14) 



and then defining the Legendre transform 



r[^,^,^n = W[M,/./l]-f drfd(7,(/J^+^V.). 

(15) 

/"[^t,^,^^] is called the effective action. In interacting 
quantum field theory F is the generator of one-particle 
irreducible Green functions [16]. It is important to 
realize that when F is defined in this conventional way, 
it corresponds to a fixed chemical potential rather than 
a fixed particle number. 
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The thermodynamical quantities might all be expressed 
in terms of the effective action. We are not going to 
present all of them but only the particle number, 



N={N) = ^iT[Nt 



-l5(H-uiN)-i 



]_J_ 



W 



fj^Jt '' 



(16) 



which is immediately seen to have the form 



A^=-^ 



]_dr 
p dfjb 



ijrijrt 



(17) 



r = In det € 



d_ 
dt' 



M 



2m 



D^'+Uxix) 



(20) 



The first two terms in Eq. (19) may be recognized as 
the classical action for the background field ^with no 
Schwinger source terms. The last term of F, which we 
have written as F , contains the effects due to the quan- 
tum fluctuations ^ ' around ^, 1/^ = 1/^+1/^'. 

We can now try to use our expression for 7^. The first 
problem we encounter is that we must obtain a more 
explicit result for F by evaluating the determinant of a 
differential operator. The most elegant method for doing 
this makes use of generalized ^-functions [18,19], and is 
motivated by analogy with the determinant of a matrix. 
If M is any Hermitian matrix, we would define det M to 
be the product of all of its eigenvalues. If m^ where 
7 = 1, . . . , n , are the eigenvalues, then 



Furthermore, in the given formulation it can be shown 
[17], that in thermal equilibrium /^[^t,^,^^] should be 
a minimum. In particular 



8^ 



= = ^ 



(18) 



To summarize this section, if we have an expression for 
the effective action F, then the effective field equations 
which determine the background field follow from 
Eq. (18). The total number of particles may be com- 
puted from Eq. (17). Although there may be easier ways 
of obtaining the basic quantities of physical interest in 
the case of free quantum field theory, the effective action 
formalism has a systematic expansion which can be used 
for interacting field theories [16]. In the next section we 
will discuss a practical way for obtaining the effective 
action. 



In det (^M) = In i^^i) = E In (^m,) . (21) 
Suppose that we define a function ^(5) by 



^(^) = EKr- 

7=1 



(22) 



By analogy with the Riemann ^-function, which is 
^r(5) = 2 ^"^ for dt{s)>\, the function defined in 

n = 1 

Eq. (22) is called a generalized ^-function. A simple 
computation shows that 

n 

aO) = n, ^'(0)= -2lnm, . (23) 

We can therefore write Eq. (21) in the form 



In det (m) = -^'(0) + ^(0) In € . 



(24) 



3. The Generalized ^-Function 

The theory described by the action Eq. (11), which 
does not involve any self-interactions, is the simplest to 
deal with because the path integral Eq. (10) defining ^ 
may be done exactly, since the integrand is just a 
gaussian. As a result, the effective action is found to be 



r=^em + 



dt da. 



2m 



iz/t 

dt 



-M- 



(19) 



where 



While this does not offer any practical advantages for 
ordinary finite dimensional matrices, it does suggest a 
possible way to define In det(^M) when M is a differ- 
ential operator: namely, set up the eigenvalue problem 
for the operator, work out the eigenvalues, and define a 
generalized ^-function as in Eq. (22). Because the 
eigenvalue spectrum of a differential operator such as 
that occurring in Eq. (20) is not bounded, in general the 
sum over all eigenvalues used to define ^(s) will diverge 
unless we restrict s to some region of the complex plane. 
This is the same as occurs for the sum used to define the 

Riemann ^-function ^r(s)= X n~' which only con- 

n = 1 

verges for ffi(s)> 1; however, the Riemann ^-function 
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may be defined by analytic continuation throughout the 
entire complex ^ -plane. (See Ref. [20] for example.) For 
the generalized ^-function we can try to define ^(0) and 
^'(0) by analytic continuation from the region of the 
complex plane where the sum over all eigenvalues con- 
verges. 

For the operator in Eq. (20) we can set up the eigen- 
value problem as follows. Let/n(x) satisfy 



2m 



D^+U,(x) 



fn(x) = EJr,(x). 



(25) 



Suppose that we define 

00 

F{s,a)= 2 {ii^j + ar (31) 

where Wj = luj/p. It was shown in Ref. [22] that 



F(s,a) = a-' + -^^ 



P' ^ e""^ 



n^)„t', «>- 



(32) 



Using this basic result, the ^-function Eq. (28) may be 
written as 



Assume that U](x) ^ so that E^ ^ 0. The eigenvalues 
denoted by E^ in Eq. (25) may be recognized as the 
energy levels for the time independent Schrodinger 
equation. We will assume that {/« (x)} forms a complete 
set of solutions to Eq. (25) which satisfies the boundary 
conditions relevant to the problem, and which are nor- 
malized by 



I 



d(r,/„*(x)/„'(x) = 5„„. 



(26) 



Because the integration over the fields ^ involved only 
those fields which were periodic in imaginary time with 
period p , the eigenfunctions of the differential operator 

r) 1 

— -fjL-:r—D^-\- Ui(x) take the form/„ (x)cxp(27TJit/l3), 
ot 2m 

and the eigenvalues of this differential operator are 



_ 27Tij 
Ajn-—^-fX-\-tn , 



(27) 



where 7 = 0, ± 1, ± 2, . . .. The generalized ^-function 
is defined to be 



as)= S S(A,J-^ (28) 



in direct analogy with Eq. (22). We will define 

T = -C(0)^aO)ln€ (29) 

again by analogy with Eq. (24). 

Before proceeding, it is helpful to show how this 
definition of the effective action makes contact with the 
standard thermodynamic results. The usual way of 
studying a system involves calculating the thermody- 
namic potential £2 defined by [21] 



)S/2=2ln 



1-e 



-j8(E„-/i) 



(30) 



^(^) = S(^«-Mr+^r(^), 



where 



^.(-)=7€^SE 



-Aj8(E„-/i) 



ns)^^, k'-^ 



(33) 
(34) 



If we expand ^t(s) about ^ = we see that ^r(O) = 
and 



n k=\ k n 



2_^-j8(E„-/i) 



(35) 



(This is noted easily from using the expansion 
\ir{s) = -I- y^^ -I- . . . valid near ^ = 0.) As T ^ we 
see that ^}(0)^0. Only the first term in Eq. (33), 
which has no^ explicit temperature dependence, will con- 
tribute to 7^ at T=0. This contribution is associated 
with the zero-point energy which arises in the path 
integral approach [23], and disappears if we normal 
order the operator H - fiN. In the ^-function method 
this normal ordering is accomphshed by taking 
as)= ^r(^). Then Eq. (29) gives 



r = -^K0) = )S/2 



(36) 



from Eqs. (30) and (35). 

To summarize this section, we have shown how the 
effective action which governs the quantum theory may 
be computed for the simple Schrodinger field theory 
described by the classical action functional Eq. (1). The 
result consists of a sum of two terms; a classical part 
involving the background field ^, and a quanturn part 
given formally by Eq. (20). The formal result for F was 
given meaning by the introduction of a generah^zed ^- 
function, and we showed how the definition of F in Eq. 
(29) was equivalent to the usual thermodynamic poten- 
tial. In the next section we will show how all of this 
formalism may be used to discuss BEC. 
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4. BEC as Symmetry Breaking 

For the free boson gas in three spatial dimensions, 
BEC can be understood as a phase transition. At the 
critical temperature characterizing the transition, the 
specific heat has a maximum, and the derivative of the 
specific heat is discontinuous. (See Refs. [24,21,25,26] 
for example.) In quantum field theory this phase transi- 
tion can be interpreted as symmetry breaking where the 
symmetry which is broken is the U(l) gauge symmetry 
associated with the change in phase of the wave func- 
tion. For charged particles coupled to electromagnetism, 
this symmetry is a local gauge symmetry; for uncharged 
particles the symmetry is a rigid, or global, symmetry. 
This was discussed in the context of relativistic field 
theory in cases where the background field was constant 
[28,27], as well as the more general case where the 
background field is not necessarily constant [29,30]. 
Nonconstant background fields are essential in cases 
where there is an applied magnetic field, or where there 
is a potential term Ui(x). 

The formalism set up in Sees. 2 and 3 is applicable 
whether the background field ^ is constant or not. The 
equations of motion for ^ and ^^ were given in 
Eq. (18). Using the result for the effective action 
obtained in Eq. (19) we find 



-^D^^U.ix)-^ 



^=0, 



(37) 



along with the complex conjugate equation. We have 
used the fact that for static potentials and electro- 
magnetic fieldsthe background field should be indepen- 
dent of time: ^ = ^ (x ). The simplicity of this result is 
also a consequence of our assumption that the theory 
does not contain any self-interactions; this assumption 
results in F containing no explicit dependence on ^ 
Symmetry breaking is associated with a nonzero value 
for ^ We may expand W in terms of the complete set 
of solutions {fn(x)} to Eq. (25): 



nx)=2cjn(x), 



(38) 



where Cn are the expansion coefficients which must be 
determined. Substitution of Eq. (38) into Eq. (37), and 
using Eq. (25), results in 



2c„(£„-m)/«(-v) = o. 



(39) 



Because the eigenfunctions obey the orthonormality 
condition Eq. (26), if we multiply both sides of Eq. (39) 
by /* (x ) and integrate over x we have 



C„(£„-M) = 0. 



(40) 



In order that the thermodynamic potential Eq. (30) 
makes sense, we must have 



fJi^ Eo 



(41) 



where Eq is the lowest energy level. In terms of the 
generalized ^-function this condition ensures that the 
effective action, or Helmholtz free energy, is real. 
It also ensures that the particle occupation numbers 
l/(e~^^^"~^^- 1) are all non- negative. It then follows 
from Eq. (40) that if ^t < £"0 ^ £"„ , the only solution to 
Eq. (40) is for all of the expansion coefficients C„ to 
vanish. In this case Eq. (38) becomes simply ^(x) = 0, 
and there is no symmetry breaking. However, if it is 
possible for the chemical potential fju to reach the critical 
value fjLc defined by 



fXc - Eq ^ 



(42) 



then Cn in Eq. (40) will be undetermined. In this case all 
of the Cn with n i= will vanish, and the background 
field is given by 



^(x) = CMx), 



(43) 



where /o(a:) is the eigenfunction corresponding to the 
ground state. If fjLc = Eq is possible to attain, the symme- 
try is broken. 

We can now make a direct link between symmetry 
breaJdng and BEC. The particle number was given in 
terms of the effective action by Eq. (17). If we use the 
result Eq. (19) for F, we may write 



N = No-^Ni 



(44) 



where 



= da, ^^ ^ 

- \r P 



and 



N,=-- 



1 dF 



P dfX 



(45) 



(46) 
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Because Co is associated with the ground state eigen- 
function in Eq. (43), it is natural to try and associate A^o 
with the number of particles in the ground state, and A^i 
with the number of particles in excited states. If we use 
Eqs. (30) and (36) we have 



iv, = E 



-j8(£„-/x)_j 



where (Oj = 2 irj/p . If we are interested in the infinite 
volume limit we may take L i , . . . , L^, to be very large, 
and replace the sums over n^, .... no with integrals. 
These integrals may be performed with the result 



(47) ^(^)= y(^)°''^^^^^,^F(^-D/2,-M) (52) 



In the next section we will see how a general criterion 
can be obtained which allows us to see if it is ever 
possible for fju to reach the critical value fjUc = Eq, and 
hence to decide if symmetry breaking can occur. For the 
remainder of the present section we will look at two 
simple applications of the results we have obtained so 
far. 



where F(s,a) was defined in Eqs. (31) and (32). Defin- 
ing ^t(s) as in Eq. (33), and making use of Eq. (36) we 
find (with T=p-') 






.npfi 



(53) 



4.1 Free Boson Gas 

The first example we will discuss is the free boson gas 
in D spatial dimensions. The special case of D = 3 is 
treated in Refs. [24,21 ,25,26] using conventional statisti- 
cal mechanical methods. For the case of general D see 
also Ref. [31]. We set A = and Ui(x) = 0. We take X 
to be a box of dimensions L\, . . . ,Ld, and will impose 
periodic boundary conditions on the field with the 
infinite box limit taken at the end. Equation (25) simpli- 
fies to 



-^V'Ux) = EJ^{x)^ 



here, and we have 



^ ^iTTfljY 



2mf__, 



Lj 



(48) 



(49) 



where n^ = 0, ± 1, ±2, . . .. The label n on En stands for 
the SQi (n I, ..., Hd) characterizing the energy levels. 
The lowest energy level is £^o = here, so that the critical 
value found for the chemical potential is fjLc = 0. We 
therefore require fjL^O. The eigenfunction/o(A:) corre- 
sponding to £^0 = is 



Mx) = v- 



(50) 



where V= Li . . . Ld is the volume of the box. (The 
factor of y'^^ comes from the normalization condition 
Eq. (26)). 

The generalized ^-function is [see Eqs. (27) and (28)] 



1 1 



'■'^^-^^i?, 



iTrriiV 



U 



(51) 



From Eq. (46) we find 






(54) 



= 1 n 



If D > 2 the sum in Eq. (54) is bounded for all ^t ^ 0. 
For large T we have A^i ~ T^^. This means that for large 
enough temperatures we can put any number of parti- 
cles into excited states. In other words, for any value of 
N, no matter how large, we can always solve N = N\ 
where N\ is given by Eq. (54) for ^t with ^K 0, provided 
that the temperature is large enough. From our discus- 
sion above this means that Co = 0, resulting in A^o = 
and ^ = so that there is no symmetry breaking. 

Now consider what happens as T decreases. As this 
happens, ^t must increase towards ^t = if we are to 
satisfy N=N\ with N\ given by Eq. (54). Eventually a 
critical temperature T^ is reached at which ^t = 0. This 
temperature is defined by 



A^ 



-^m 



UD/2) 



(55) 



where ^r is the Riemann ^-function. If p = N/V is the 
density of particles, we have 



r.= 



277 



r 1^'^ 



For r < Tc it is not possible for ^t to decrease beyond 
^t = 0, so it remains frozen at this critical value. From 
Eqs. (54) and (55) we have 



n.-n[L)' 



(57) 
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It is therefore not possible to accommodate all of the 
particles in the excited states. Because A^o = N-N\, we 
have 



A^o=A^ 



(58) 



as the number of particles in the ground state. Using Eq. 
(45) we find 



(59) 



4.2 BEC in a Magnetic Field 

We will now consider the case where there is a 
constant applied magnetic field. We will allow the 
spatial dimension D to be general, but will only consider 
the case of a magnetic field with a single nonzero com- 
ponent. (The general case is more complicated and is 
discussed in Ref. [17].) If we choose the magnetic field 
to point in the z -direction, a suitable gauge choice for 
the vector potential is 



A = (-53;,0,...,0) 



(64) 



and from Eqs. (43) and (50) we have for T < Tc 



where B is the strength of the magnetic field. Eq. (25) 
reads 



^ = p' 



'-'^ 



(60) 



For the special case Z) = 3, the results Eqs. (56) - (58) 
reproduce the standard expressions. The generalization 
to arbitrary D was given by May [32]. 

If Z) ^ 2 the situation is entirely different than the one 
we have just described. If we try to let fx ^ in 
Eq. (54), it is observed that the sum is not bounded. This 
means that we may put any number of particles into the 
excited states. Equivalently, for Z) ^ 2 we can always 
solve N = N\ with A^i given in Eq. (54) for ^ with ^t < 
for any temperature. It is not possible for ^ to reach the 
critical value ^t = for any T> with a finite number 
density of particles. This is easy to see when D = 2 
because the sum in Eq. (54) may be easily performed to 
give 

iVi=-v(f|)ln[l-e^^]. (61) 



2m \dx J 



1 f d' 

' 2m f^^ b(x'f 



fn(x) = EJ„(x). 



(65) 



This equation is equivalent to that for a simple harmonic 
oscillator, and the energy levels are easily found to be 

[34] 



eB 1 ^ /27772A2 



eB 1 ^ /lirnA 



(66) 



where 7 = 0, 1, 2, ... ; rii =0, ±1, ±2, . . . ; and we 
have again imposed periodic boundary conditions. The 
eigenvalues Eq. (66) are degenerate with degeneracy 
eBLiL2/(27T). The smallest energy eigenvalue is seen to 
be 



2m 



(67) 



Setting A^i =N= pV, and solving Eq. (61) for ^ gives 
/^ = rin[l-e-^]. (62) 



and the critical value of the chemical potential is 



Mc = 



eB^ 

2m ' 



(68) 



For small T this results in 



The generalized ^-function defined in Eq. (34) is 



Jlirp 
fJi— OTrnT . 



(63) 



eB T 

Ct{s)= 1Z—L1L2- 



27r ' 'r(s) 



There is no BEC, in the same sense as BEC for the Z) = 3 
gas, when D = 2. This agrees with the analysis of 
Refs. [32,33]. For Z) = 1 it is not possible to perform the 
sum in Eq. (54) in terms of simple functions for ^ < 0; 
however, since the sum is not finite for ^t = 0, BEC is not 
possible in this case. 



s s i s£^. (69) 



; = n3 = -oo n£, = -oo k=l 



Taking the large box limit, and replacing the sums over 
^3, . . . , ^D with integrals (which just involves a product 
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of gaussians), and noting that the sum overy is just a 
geometric series, results in 



^t(s)= V 



eB\ (mT 



477/ \ 277, 



E 



-npfx 



^ (^) n=\ ^D/2-s „:„!_ fneB 



n^'^-^sinh 



2mT^ 



(70) 



State, then it may still occur even in the absence of a 
phase transition. It is important to keep the definition 
which is chosen for BEC firmly in mind. We will return 
to this matter in Sec. 6. Finally, it can be shown that even 
though there is no phase transition, the charged Bose 
gas still exhibits a Meissner-Ochsenfeld effect [35]. 
Although we have shown that there is no phase transi- 
tion if Z) ^ 4, we have not shown that there is one for 
D > 4. To do this we must show that A^i remains 
bounded as fji ^> fjic. If we use the inequality 



From Eq. (36) we find 



r = -a{0) = -v 



eB\ (mT 



1 



477/ \ 277, 



« = 1 ,^D/2 



n^^sinh 



neB 
2mT/ 



(71) 



Using this result in Eq. (46) leads to 



A^i 



-(i)(i: 



-n^iix-ix^) 



neB 

l-e~mr 



(72) 



with [jLc given by Eq. (68). 

We can now analyze whether or not BEC occurs in 
the same way as for the free boson gas in the previous 
example. First, if we use the inequality (l-e""")"^ > 1, 
valid for all x > 0, we see that 






(73) 



From our earlier discussion we know that BEC is only 
possible if A^^i remains bounded as /x ^ fic. Because the 
sum in Eq. (73) is not bounded as fi ^ fic for D/2- 
1^1, this inequality shows that BEC will not occur for 
D ^ 4. This includes the physically interesting case of 
Z) = 3, as shown originally by Schafroth [35]. The 
absence of BEC for Z) ^ 4 was given originally by May 
[36]. 

Of course one must be precise about what is meant by 
the absence of BEC here. What we have shown is that if 
BEC is interpreted to be synonymous with symmetry 
breaking and a phase transition in the same way as BEC 
occurs for the free boson gas in three dimensions, then 
it does not occur. On the other hand, if one interprets 
BEC as a sudden build-up of particles in the ground 



(neB\ / eg \ 

1-e mr 1 ^ I 1-e mrl 



valid for all n ^ 1, then we see that 



(74) 






■n^(jLt-jLtc) 



(75) 



For Z) ^ 5, the sum on the RHS of this inequality 
remains bounded as /x ^ fjic. As in the example dis- 
cussed in Sec. 4.1, this means that it is not possible to 
place an arbitrary number of particles in excited states 
if r < Tc where T^ is the solution to 



A^ 



^(mmr i'-i^-'^-j 



(76) 



Unlike the case of the free boson gas, it is not possible 
to obtain an explicit expression for Tc in terms of the 
particle number, although approximate expressions can 
be obtained for strong and weak fields [17]. 

To summarize this section, we have shown the con- 
nection between BEC and symmetry breaking. The gen- 
eral results were illustrated with two examples. The first 
was the familiar case of a gas of free bosons. The second 
was a gas of charged bosons in a constant magnetic 
field. In both cases we saw how BEC interpreted as a 
phase transition corresponding to a breaking of the t/(l) 
gauge symmetry could occur in some cases but was 
inhibited in others. The crucial deciding factor was the 
number of spatial dimensions. In the next section we 
will discuss a general criterion for deciding if BEC as 
symmetry breaking can occur. 

5. General Criterion for BEC as 
Symmetry Breaking 

In the last section we saw how the free boson gas 
did not undergo BEC, at least in the sense of a phase 
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transition, if the spatial dimension Z) ^ 2. If a constant 
magnetic field was applied to a charged gas of bosons, 
no BEC occurs if Z) ^ 4. In this section we will discuss 
the underlying features present in these two examples 
which allows a unified treatment of these two cases, and 
in addition generalizes the analysis to a wide class of 
systems. We follow Refs. [37,38]. 

The key feature present in both examples given in the 
last section is that the energy levels may contain a dis- 
crete part as well as a continuous part. For the free Bose 
gas in D dimensions we had Eq. (49). In the infinite box 
limit /2i, . . . , /2£, could be treated as continuous. There- 
fore £'„was labelled by D continuous quantum numbers. 
BEC was found to occur for Z) ^ 3. For the charged 
Bose gas in a constant magnetic field the energy levels 
were given by Eq. (66). This time, again in the large box 
limit with fi3, . . . , n^, treated as continuous, the energy 
levels involved (D-2) continuous quantum numbers, and 
one discrete label corresponding to the degenerate 
Landau levels. This time BEC occurs if Z) ^ 5, which 
can be suggestively written involving the number of 
continuous quantum numbers as (D-2) ^ 3. The 
feature common to both examples given in Sec. 4 is that 
the dimension of the space associated with the continu- 
ous labels had to be at least 3 for BEC to occur in the 
sense of a phase transition and symmetry breaking. 

Suppose that we consider any system for which the 
energy levels can be expressed as the sum of a discrete 
part which we will denote by Ep, and a continuous part 
which we deal with by box normalization with the 
infinite box limit taken at the end as in the examples 
presented in the last section. We will assume that the box 
has dimension q . We may write 



rLfi — rLr, + 



2m 



277/1. 



(77) 



where Li , . . . , L^ are the sides of the box. Here/7 is just 
a set of labels for the discrete part of the spectrum. With 
the large box limit taken, the labels /2i, . . . , n^ may be 
regarded as continuous. The generalized ^-function (34) 
reads 



^t(s) = 



fis) 



1 fd^-i 

p J k=i 



,-kP(E„-n) 






(78) 



after the integration over the continuous part of the 
energy spectrum has been performed. 

We found that BEC with the associated phase transi- 
tion only occurs if it is possible for ^t to reach the critical 



value |uic determined by the lowest energy level. (See the 
discussion around Eq. (42).) In the case of Eq. (77) we 
have fjic = Eq = Eq. Because the lowest mode is playing 
such a crucial role, it is advantageous to separate it off 
and define 



^^(^)=^r(*)+ ^r\s) 



where 



^r(*): 



Vgdo T""--' 



-t^(Mc-M) 



(477-)'^ r{.s) t, V'"^- 



(79) 



(80) 



represents the contribution of the lowest mode to the 
^-function (with do the degeneracy), and ^t^^^ {s) is 
given by Eq. (78) but with the sum over/7 restricted to 
the nonzero modes only. We have ^t (0) = as before, 
and 



^f '(0) = 



Vado 



(471) 



q/2 



r^/^2 



-kp(tic-t^) 



^l+,/2 ' 



(81) 



^r"(0): 



(477-) 



q/2 



r'^S E^t^t::;^- (82) 



p#0 k= 1 



^l+q/2 



These expressions may now be used to find the effective 
action F or the thermodynamic potential O given in 
(36). 

The advantage of separating off the lowest mode as 
we have done is that the argument of the exponential in 
Eq. (82) will remain negative even if ^t = fjLc. The conse- 
quence of this is that the sums in Eq. (82) will converge. 
This remains true even if we differentiate Eq. (82) with 
respect to fjL to find the contribution of the^excited states 
to the particle number given in terms of F by Eq. (46). 
Therefore, whether or not a phase transition occurs rep- 
resenting BEC is determined by the behaviour of ^^^\s). 
If we drop terms which remain finite as ^t ^ ^tc we 
have 



_d_ 
dfX ' 



N(fJL^fjL.)^T^^P'(0) 



J- \q/2 



Mo(£) 



k=\ 



-k(S(iXc-ix') 



k'"^ 



(83) 



If the sum in Eq. (83) remains finite as ^t ^ ^tc, then 
the number of particles which can exist in excited states 
is bounded, and we get BEC and a phase transition with 
symmetry breaking. Clearly this can occur only for 
^ > 2. If i^f ^ 2, we can conclude that BEC does not 
occur, at least in the sense of a phase transition. 
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We can obtain much more detailed information on 
exactly how A^ diverges as ^t — > ^t^ ■ For q = the sum in 
Eq. (83) is just a geometric series, and for ^ = 2 the sum 
is just the expansion of the natural logarithm. For q= I 
the sum may be evaluated as described by Robinson 
[39]. We therefore find 



N(fJL -^ fJLc) 






(for q = 0); 



(84) 



N(fi ^ Mc) - 2 ^^'"^^ (Mc-m)"^'' (for q = 1); 



(85) 



N(fJL -^ fJLc) 



1 / T 
-TV2do\n{ 



(for q = 2). 



(86) 



Only that part of A^ which diverges as fx ^ fXc has been 
shown. The ground state when there is symmetry break- 
ing can be determined in the way described in Sec. 4. 
In this section we have shown how the occurrence of 
BEC is linked to the number q which characterizes the 
continuous part of the energy spectrum. For BEC to 
occur we require q ^3. For the free boson gas in D 
dimensions, since q = D, this recovers our earlier result 
in Sec. 4.1. For the charged boson gas in a constant 
magnetic field, we have q = D-2, again recovering our 
earlier result in Sec. 4.2. The physical meaning of q is 
that it is the number of spatial dimensions for which the 
particles are effectively free to move. For the boson gas 
in a magnetic field in the z direction, classically the 
motion in the x-y plane is restricted to circular orbits, 
and it is only in the z -direction that the motion is free. 
A variety of other examples, often established by long 
and detailed calculations all emerge from this relatively 
simple analysis [38]. A similar approach may be used 
for relativistic gases [37,38]. 

6. BEC in Harmonic Oscillator 
Confining Potentials 

We showed in the last section how BEC interpreted as 
symmetry breaking could occur only if the continuous 
part of the energy spectrum involved at least three con- 
tinuous labels (q ^ 3). A special consequence of this is 
that if the energy spectrum is entirely discrete, corre- 
sponding to ^ = 0, then BEC as symmetry breaking will 
never occur. In other words, for a system characterized 
by a discrete set of energy levels, if BEC does occur it 
cannot be interpreted as a phase transition analogously 
to the free boson gas in three dimensions. However, as 
we have discussed recently [10], it is still possible to 
have BEC in the sense that there is a critical temperature 
characterizing the system at which the number of parti- 



cles in the ground state has a sudden and dramatic 
increase (see also Ref. [40]). This is borne out in the 
experiments on alkali gases at microkelvin temperatures 
[7-9]. 

The simple model we will discuss here is a system of 
uncharged spin-0 bosons in a harmonic oscillator con- 
fining potential. This is of phenomenological interest 
since it represents a model for the magnetic traps used 
in experiments [7-9]. We will use the action Eq. (1) with 
A = 0, and U\(x) the harmonic oscillator potential 



Ui(x) = -m((oh^-^ (oiy^-\- ojW) - 



(87) 



(We will only consider 3 spatial dimensions here 
although the analysis may be generalized in an obvious 
way to any number of dimensions.) The energy levels 
are simply 



E^ = (ni -I- x)^wi -I- (^2 + 2)^^2 + (n3 + 2)^<^3 , (88) 



where n = {n\,n2, n^,) with n, = 0, 1, 2, . . . and we have 
reinstated the explicit /i dependence. The lowest energy 
level, which determines the critical value of the chemi- 
cal potential as in Eq. (42), is 



fji^ = Eo= jfi{o}x -I- W2 + o^z) ■ (89) 



The total number of particles is given by the usual Bose 
distribution function 



A^ 



= i: i i |e^^--->-ir 

n 1 =0 n2 = n 3 = '^ J 



(90) 



Although it is possible to proceed with the anisotropic 
case, our results are most easily illustrated for the 
isotropic harmonic oscillator with (Oi = 0)2= 0)3 = (o. 
(The general anisotropic case is dealt with in Ref. [1 1].) 
For the isotropic oscillator, the triple sum in Eq. (90) 
becomes the simpler result 



A^=S|(/+l)(/ + 2){e^^"^^^-l} 



(91) 



where we have defined the dimensionless variables x 
and € by 



kT 



fX=^0) (--€). 



(92) 



(93) 
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Since ^^ = l^(o, a phase transition with the associated 
symmetry breaking occurs only if it is possible for fju to 
reach the value 6=0. Our general analysis of the last 
section shows that this can never happen; however, for 
this particular example we can show this another way. If 
we expand 



,(/ + e)x 



1 =2e-"^^"^ 



(94) 



then the sum over / in Eq. (91) is easily done with the 
result 



of integration in an appropriate way obtain at least an 
asymptotic expansion for some appropriate range of the 
parameters. The details of this procedure are described 
in Ref. [11]. We will illustrate the general technique for 
the number of particles as given in Eqs. (91) or (95). 

The aim is to obtain an asymptotic expansion of A^ 
valid for small x and small e. We can do this by making 
use of the Mellin-Barnes representation for the exponen- 
tial: 



-^ = y^ dar(a)v- 



(97) 



A^= 2e""'"(l-e-"T'- (95) 

n= 1 

For X > we have (l-e"""")"^ > 1, so that 



A^> 2e"""=(e"^-l)- 



(96) 



As we let 6 — > 0, which should be the signature for BEC, 
it is seen that A^ increases without bound. This is quite 
different from the free boson gas in three dimensions 
where the number of particles is bounded. This means 
that regardless of the temperature, it is always possible 
to solve Eq. (95), or equivalently Eq. (91), for e given 
any finite N, with 6>0. It is never possible for this 
system to attain the limit 6 = for any T > and any 
finite number of particles. The fact that a Rose gas in a 
harmonic oscillator potential does not condense in the 
same way as a free boson gas was noted originally in 
Ref. [41]. 

This points out a fundamental difference between our 
analysis and other treatments [42,40,43]. We have 
treated the energy spectrum for the harmonic oscillator 
as discrete with the ensuing sums. The physically inter- 
esting case occurs when (o/lv— 100 Hz, and T -- fxK. 
In this case, x as defined in Eq. (92) is small. A plausible 
approach is to argue that for x« 1, it is justified to 
replace a sum such as Eq. (95) with an integral. This is 
tantamount to regarding the energy levels as continuous 
rather than discrete, and the analysis of Sec. 5 shows that 
the underlying physics is crucially dependent on the 
number q of continuous dimensions in the energy spec- 
trum. Any approximation which effectively changes 
q from to 3 is therefore suspect. 

If the correct behaviour for small x is desired, the only 
safe approach is to deal with the exact sum Eq. (95). The 
result in Eq. (95) does not converge very rapidly for 
small X, nor does it display in any transparent way the 
behaviour at small x. However, it is possible to convert 
Eq. (95) into a contour integral (which is an exact result, 
not an approximation), and by deforming the contour 



valid for 3?(i^) > and c E ^ with c> 0. Equation (97) 
is easily proven by closing the contour in the left hand of 
the complex plane, enclosing the simple poles of r(a) 
at q; = -fz , fz = 0, 1 , 2, . . . with residue ^^-p . The 
residue theorem immediately gives the Maclaurin series 
for e~^ From Eqs. (91) and (94) we have 



N= E E^(/+l)(/ + 2)e-"^^^^^- 






(98) 



The order of summation and integration may be inter- 
changed provided that we deform the contour of integra- 
tion first so that c > 3. The sum over n can now be done 
in terms of the Riemann ^-function, and the sum over / 
in terms of the generalized, or Hurwitz, ^-function 
^H{s,a) defined by Ref. [20] 



^fj(s,a)= ^ (n-\-a) 



(99) 



for 9i(s) > 1 and 0<a ^ 1. (The Riemann ^-function is 
the special case ^r(s) = ^h(s, 1).) The sums in the orig- 
inal expression for the particle number have now been 
performed exactly with no approximations, and we 
have obtained an integral representation for the particle 
number: 



A^ 






+ (3-26)^^(a-l, €) + (1-6) (2-6)^^(a, e)} . 

(100) 
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Closing the contour in the left half of the complex plane 
and use of the resldure theorem now leads to an asymp- 
totic expansion of A^ in powers of x : 



j,^^Jl^U'.-eV-^.^.0M. (101) 



It is possible to extend the range of validity of the expan- 
sion to larger values of e. This can be done by treating 
the first (/ = 0) term in Eq. (91) separately and applying 
the contour integral procedure we have just described to 
the remaining terms. The result is 



N'- 



1 ^/.(3) /3 \U2) ^/in 



(102) 



which is more accurate than Eq. (101) for larger €. 
Similar expansions may be obtained for other thermody- 
namic functions [10,1 1]. For the internal energy we find 



U 3^«(4) (9 , U«(3) 13 ^« (2) 
■7— = — +k-3e — — + ■ 



4x' 



2 ex \x X 



For the specific heat we find 



(103) 



^,, 12^/^(4) 9^r(3) 2^r(2) 126^^,(3) 
C/k = + + ^Ji^-^ 



9€'^A3f \S€'^A2)^r(3) 



, 9.Y.gK.(3)^ ^(,^^^£)^ 



(104) 



These analytic results can be shown to be in excellent 
agreement with numerical results obtained from an eval- 
uation of the exact sums in the range in which they are 
valid ( small x and small e). 

In a similar way we are able to treat the anisotropic 
oscillator potential. Introducing the geometric mean 
n = (wi (02(02y^^ of the oscillator frequencies, we found 
for the particle number 



AT^- 



1 



,(;5;ee/3)(wi+W2+W3) _ 



-A0^'^'^''{m)^'^'^^ 



with 



2 \(Oi( 



1 1 

■ + + ■ 



0*2 (Oi(03 0)2(03 



(105) 
(106) 



Analogous results for U and C are given in [11]. 

Because there is no sharp phase transition, identi- 
fying a critical temperature is problematical. One 
approach which has been used in finite volume systems 
[44,45] is to calculate the maximum of the specific heat, 
and identify the temperature at which the maximum 
occurs with the BEC temperature. It is very difficult to 
obtain reliable analytic expressions for the specific heat 
in the region of the maximum for this model. However 
it is possible to compute the specific heat numerically, 
and the result is shown in Fig. 1. 

We have chosen ft;/27r= 416 Hz to be the geometric 
mean of the frequencies in the sodium experiment [9]. 
The number of particles A^ = 5 X 10^ is also taken from 
the sodium experiment. The maximum in the specific 
heat occurs forx — 0.0136, corresponding to a tempera- 
ture of r — lAl fjL K. This is in remarkably good agree- 
ment with the temperature of 2 ^tK quoted in Ref. [9]. 

Another difference between our results and earlier 
work [42,43] is that we find the specific heat to be 
smooth and continuous at its maximum. A closeup of 
the specific heat in a neighbourhood of its maximum is 
shown in Fig. 2. The discontinuous behaviour found in 
Refs. [42,43] is due to approximating sums with inte- 
grals, which as stated earlier is not a reliable approxima- 
tion. 

That the specific heat is continuous has also been 
realized by the authors of Ref. [40]. The question of 
whether or not one can actually distinguish in an exper- 
iment a drop as seen in Fig. 2 from a genuine disconti- 
nuity has also been addressed there. They improved the 
analysis of [42,43] by taking into account a suitable 
density of states. In their procedure a parameter depend- 
ing on the oscillator frequencies had to be determined 
numerically. Our above described procedure gives the 
complete analytical dependence, see Eq. (106). 

Another possibility for obtaining an estimate of the 
critical temperature for BEC consists of examining the 
population of the ground state. The first term on the 
RHS of Eq. (102) may be observed as the number of 
particles in the ground state. (Put ^ = in Eq. (91).) The 
remaining terms then give the number of particles in 
excited states. We could define the critical temperature 
to be the temperature at which a specified fraction of the 
total number of particles are in the ground state. If we 
define 



where 



7Vo=(e^^-l)-', 



(107) 
(108) 



this fixes € , and hence the chemical potential, in terms 
of X and N by 



€X 



='"('^i^) 



(109) 
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Fig. 1. The specific heat in units of/: as a function of x = fi(i)l{kT). The total number of particles 
isA^ = 5XlO^ and w/27r = 416Hz. 
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Fig. 2. The specific heat in units of A: as a function of lOOx where x = fi(t)l{kT). 



For a large number of particles, e would be expected to 
be small. If we use Eq. (102) we also have 



{\-f)N^^,0)x-' + ^^,{2)x-' 



(110) 



which gives a cubic equation for x . This may be solved 
in a straightforward way. (A more accurate result may be 
obtained by including more terms in the asymptotic 



expansion [Eq. (102)] beyond those indicated [11]). 
With iV= 5X 10' we find x ^ 0.0136 for /= 1/100; 
x=- 0.014 for /=1/10; x =- 0.017 for /=0.5; and 
X — 0.0373 for/= 0.95. For the purpose of comparison, 
we have computed the ground state population numeri- 
cally. The result is shown in Fig. 3. The approximate 
result for the number is in very good agreement with the 
exact result. 
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500 




Fig. 3. The number of particles in the ground state as a function of x. 



We have worked out the maximum in the specific heat 
for the two cases A^= 2 X 10^ and A^= 2 X 10^ and find 
that it corresponds to x=^ 0.0185 and x =^ 0.0408, 
respectively. Using the oscillator frequencies given in 
Ref [8] we find the temperature at which the specific 
heat has a maximum is T — 380 nK, in good agreement 
with the range of 100 nK to 400 nK for the experiment. 
For the case of rubidium with N=2X\^'^ we find 
r — 7 1 nK if we use the oscillator frequencies of 
Ref. [7] . If we use the frequencies given in Ref. [46] for 
the strong trap we find T— 124 nK, again in close 
agreement with the experiment. We can also compare 
the results for the maximum in the specific heat to the 
fraction in the ground state found from Eq. (110). For 
the case A^= 2X 10^ we find x ^ 0.0185 for/= 0.01; 
X- 0.0191 for /= 0.1; x - 0.0233 for /= 0.5; and 
X =- 0.051 for/= 0.95. If we take iV= 2X 10^ we find 
xc- 0.0404 for /=0.01; x - 0.0417 for /=0.1; 
X =- 0.051 for/= 0.5; and x =- 0.114 for/= 0.95. In all 
three cases, the maximum in the specific heat occurs 
when only about 1 % of the particles are in the ground 
state. The specific heat maximum is therefore a good 
indicator of the onset of BEC. 

For the case that BEC in the sense of a phase transi- 
tion occurs, the critical temperature is the temperature at 
which the ground state starts to become populated. (See 
Eq. (58) for the free boson gas.) Because we have seen 
that the specific heat maximum also corresponds to the 
point at which the ground state population begins to 
grow, we believe that this gives a good and reliable 
indicator for the onset of BEC in neutral atoms trapped 
by a confining potential. 
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